Dynamics of an anisotropic Haldane antiferromagnet in strong magnetic field 
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We report the results of elastic and inelastic neutron scattering experiments on the Haldane- 
gap quantum antiferromagnet Ni(C5Di4N2)2N3(PFe) performed at mK temperatures in a wide 
range of magnetic field applied parallel to the S — 1 spin chains. Even though this geometry 
is closest to an ideal axially symmetric configuration, the Haldane gap closes at the critical field 
H c ~ AT, but reopens again at higher fields. The field dependence of the two lowest magnon modes 
is experimentally studied and the results are compared with the predictions of several theoretical 
models. We conclude that of several existing theories, only the recently proposed model [Zheludev 
et al., cond-mat/0301424 ] is able to reproduce all the features observed experimentally for different 
field orientations. 

PACS numbers: 75.50.Ee,75.10.Jm,75.40.Gb 



INTRODUCTION 



The problem of field-induced magnon condensation 
in gapped quantum antiferromagnets is currently re- 
ceiving a great deal of attention from experimental- 
ists. Particularly important results were obtained in re- 
cent neutron scattering and ESR measurements on S = 
1 Haldane-gapi^ compounds Ni(C 2 H 8 N2)2N0 2 (C104) 
(NENP), 3 . Ni(C 5 Di4N 2 )2N 3 (PF 6 ) (NDMAP) 4 i 5 i 6 i 7 : 8 : 9 
and Ni(C 5 Hi4N 2 )2N 3 (C104) (NDMAZ), 10 as well as the 
S = i-dimer system TlCuClg. 11 ' 12 : 13 The effect of mag- 
netic field is to drive the gap A in such systems to zero 
by virtue of Zccman effect, thus promoting a quantum 
phase transition to a new magnetized state at some criti- 
cal value of the applied field H c w (A/g/xs), 14 Additional 
magnetic anisotropy effects usually lead to more complex 
behavior and richer phase diagram. Anisotropy is negli- 
gible in many 5* = ^-based materials such as TICUCI3, 
where no single-ion terms are possible. In contrast, for 
S = 1 compounds such as NENP, NDMAP and NDMAZ, 
terms of type DS^ are quite strong and the corresponding 
zero-field anisotropy splitting of the excitation triplet is 
comparable in magnitude to the Haldane gap itself. Un- 
der these circumstances, the physics is expected to de- 
pend strongly on the direction of the applied magnetic 
field with respect to the anisotropy axes^ 

For purely technical reasons, in quasi-lD materials it 
is much easier to perform inelastic neutron scattering ex- 
periments in high magnetic fields applied perpendicular 
to the spin chains. On most instruments the scatter- 
ing plane is horizontal and the wave vector resolution 



along the vertical axis is deliberately coarsened to pro- 
vide an intensity gain. To optimize wave vector resolu- 
tion along the spin chains in the sample, and to allow 
the momentum transfer qii w tt in that direction, the 
chain axis has to be mounted in the horizontal plane. 
The typical construction of superconducting magnets is 
such that the field is along the vertical direction, and 
is thus applied perpendicular to the spin chains in the 
sample. In the Haldane-gap materials NENP, NDMAP 
and NDMAZ the anisotropy easy plane is roughly per- 
pendicular the chains. As a result, most of the previous 
neutron measurements were performed for in-plane mag- 
netic fields, i.e., in the Axially Asymmetric (AA) geom- 
etry. In the AA case the transition at H c is expected to 
be of Ising type^^ and even an isolated chain acquires 
antiferromagnetic long-range order in the ground state at 
H > H c . The excitations in the magnetized state are a 
triplet of massive "breathers" (soliton-antisoliton bound 
states) Recent neutron scattering studies of NDMAP 
in the AA geometry provided a solid confirmation of these 
theoretical predictions. 9 

For an isolated Haldane spin chain magnetized in 
the Axially Symmetric (AS) geometry (with a magnetic 
field applied parallel to the anisotropy axis and the 
anisotropy being of a purely easy-plane type) theory pre- 
dicts a totally different, disordered and quantum-critical 
ground state whose low-energy physics can be described 
as the Tomonaga-Luttinger spin liquid»iS* 1 ii2*i 9 " The low- 
energy excitation spectrum contains no sharp modes and 
is instead a continuum of states, much as for S = o 
spin chains. 20 Higher modes which would have quasipar- 
ticle character in the AA geometry (two upper members 
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of the Zeeman-split triplet) also develop into continua 
and exhibit only edge-type singularities in the AS case^i 
Moreover, one expects incommensurate correlations at 
a field-dependent Fermi wave vector that characterizes 
the Fermi sea of magnons "condensed" at H > H c . 22 
In a real material this idealized picture may be com- 
plicated by several factors. First of all, there can be 
additional anisotropy terms, such as in-plane single-ion 
anisotropy or Dzyaloshinskii-Moriya interactions, which 
explicitly break the axial symmetry and favor the AA 
physics; a similar effect could be expected if the direction 
of the applied field deviates slightly from the symmetry 
axis. Residual three-dimensional inter-chain couplings 
lead to a spontaneous breaking of the axial symmetry 
which is equivalent to the Bose-Einstein condensation 
of magnons >ii*i2ii2i2^ Due to the critical nature of the 
ground state of an isolated chain, all these effects are 
relevant and will have a significant impact on the spin 
correlations, no matter how small they may be. The 
key to understanding the high-field behavior in real ma- 
terials is a combined approach involving high resolution 
neutron scattering experiments and a consistent theoret- 
ical treatment. Due to geometrical constraints described 
in the previous paragraph, high-field neutron measure- 
ments in the AS configuration are technically challenging 
and require the use of specialized neutron instruments or 
magnet environments. 

The purpose of the present paper is twofold: on 
one hand, to study experimentally the response of the 
NDMAP system in the wide range of applied fields and 
at very low temperature, in the setup as close to the AS 
geometry as possible, and to test further the phenomeno- 
logical field theory of the high-field phase which was re- 
cently proposed^ and successfully applied to the descrip- 
tion of experiments on NDMAP in the AA geometryti 
On the other hand, we also give a more detailed account 
of the theory which was only briefly outlined in Ref. 0, 
and generalize it to include the effect of interchain inter- 
actions. Finally, we perform a systematic quantitative 
comparison between our experimental findings and theo- 
retical results. 

The paper is organized as follows: in Sect.[n]we present 
the details of experimental setup, and Sect. Hill reports 
the results of elastic and inelastic neutron scattering mea- 
surements. In Sect. IIVI we present the general effective 
field theory of anisotropic gapped quasi-one-dimensional 
spin system in strong magnetic field and perform a sys- 
tematic quantitative comparison between our experimen- 
tal findings and theoretical results. Finally, Sect. Ivlcon- 
tains the summary and concluding remarks. 



II. EXPERIMENTAL SETUP 

A. Structural considerations and energy scales 

The crystal structure of NDMAP is schematically 
shown in Fig. 1 of Ref. 0- The AF spin chains are 
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FIG. 1: The H — T phase diagram of NDMAP as deduced 
from specific heat measurements (open symbols, Ref. [a) and 
neutron scattering experiments (solid symbols, Refs. ITU and 
the present study). The shaded area is the high-field phase 
investigated in the present work. 



composed of octahedrally-coordinated S — 1 Ni 2+ ions 
bridged by azido-groups. The chains run along the c 
axis of the orthorhombic structure (space group Pnmn, 
a = 18.046 A, b = 8.705 A, and c = 6.139 A). Previous 
zero-field neutron studies provided reliable estimates for 
the relevant magnetic energy scales in the system. The 
in-chain exchange constant is J = 2.6 meV. Exchange 
coupling along the crystallographic b axis is considerably 
weaker, Jb/J ~ 10~ 3 , and that along a is weaker still, to 
the point of being undetectable: \J a /J\ < 10~ 4 . Mag- 
netic anisotropy in NDMAP is predominantly of single- 
ion easy-plane type with D / J sa 0.25. In addition, there 
is a weak in-plane anisotropy term of type E(S 2 — S 2 ). 
As a result of these anisotropy effects, at zero field the 
degeneracy of Haldane triplet is fully lifted and the gap 
energies are Ai = 0.42(3) meV, A 2 = 0.52(6) meV, and 
A3 = 1.9(1) meV. Correspondingly, the critical fields are 
strongly dependent on field orientation and the H — T 
phase diagram, visualized in Fig. ^ is highly anisotropic. 

The local anisotropy axes in NDMAP are determined 
by the geometry of the corresponding Ni 2+ coordination 
octahedra and do not exactly coincide with crystallo- 
graphic directions. Instead, the main axes of the NiNg oc- 
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FIG. 2: Main elements of the NDMAP crystal structure in 
projection onto the (ac) crystallographic plane. Only the N 
and Ni sites are shown. Magnetic anisotropy associated with 
the Ni 2+ magnetic ions is determined by the local symmetry 
of NiN6 coordination tetrahedra that are tilted in the (ac) 
plane by a = 16° relative to the c axis. The tilt direction is 
opposite in adjacent chains. 

tahedra are in the (a, c) crystallographic plane, but tilted 
by a w 16° relative to the c axis. As illustrated in Fig. El 
within each chain, the tilts are in the same direction for 
all Ni 2+ -sites, so there is no intrinsic alternation in the 
chains, as is the case in many related compounds such as 
NENP. However, within the crystal structure there are 
two types of chains related by symmetry, and the corre- 
sponding tilt directions are opposite. This circumstance 
has a very important consequence for the present study. 
It implies that for NDMAP one can not apply the field 
along the anisotropy axis of all the magnetic ions in the 
sample. The closest one can come to this idealized AS 
scenario is by applying a field along the axis of bulk mag- 
netic anisotropy, i.e., along the (0,0, 1) direction. In this 
case the field will form a small angle of ±a with the local 
anisotropy axes for all spin chains. 

B. Experimental procedures 

Single-crystal neutron scattering experiments in mag- 
netic fields applied parallel to the spin chains were car- 
ried out using two different setups. A series of diffrac- 
tion measurements was carried out using a vertical-field 
6 Tesla cryomagnet installed on the rather unique D23 
lifting counter diffractometer at ILL. On this instrument 
the data collection is not restricted to a given scattering 
plane. The sample is mounted with the chain axis ver- 
tical (parallel to the field) and the detector is lifted out 
of the horizontal plane to allow a momentum transfers 
<7H w 7r in that direction. Neutrons of a fixed-incident 
energy Ei = 14.7 mev were provided by a curved ther- 
mal neutron guide with supermirror coating, a Pyrolitic 



Graphite (PG) monochromator and a PG filter. In some 
cases horizontal and/or vertical collimators with a 40' 
FWHM beam acceptance were inserted in front of the 
detector. In the diffraction study we employed a w 0.3 g 
fully deuterated single crystal sample. 

Inelastic measurements were performed using the con- 
ventional 3-axis cold neutron spectrometer FLEX in- 
stalled at HMI. An assembly of 3 deuterated single crys- 
tals with total mass about 1 g were mounted with chain 
axis in the horizontal (scattering) plane of the instru- 
ment. A magnetic field was applied along that direction 
by a horizontal-field cryomagnet. The data were collected 
with the final neutron energy fixed at 5 meV. A Be fil- 
ter was used after the sample to eliminate higher-order 
beam contamination. Beam divergencies were defined by 
the critical angle of the cold-neutron guide, and by char- 
acteristics of the PG analyzer and monochromator used. 
No additional devices were used to collimate the neu- 
tron beams. In both the D23 and FLEX experiments the 
sample environment was a 3 He- 4 He dilution refrigerator. 



III. RESULTS OF NEUTRON SCATTERING 
MEASUREMENTS 

A. Magnetic long range order 

Applying a field H > H c w 4.0 T parallel to the c 
axis at T — 25 mK leads to antiferromagnetic long-range 
ordering in NDMAP. This was deduced from the appear- 
ance of new Bragg reflections at the 3D AF zone-centers 
(h, k, I) + (0, —0.5, 0.5), where h, k and I are integer. The 
measured field dependence of the (0,1.5,0.5) background- 
subtracted peak intensity is plotted in Fig. [3] In this data 
set the experimental error bars are too large to allow an 
accurate determination of Tjy and the order-parameter 
critical index. 

Fig. 21 shows h, k and I scans across the (0,1.5,0.5) 
peak measured at H = 5.8 T. These data were taken 
with a horizontal collimator installed in front of the de- 
tector, to improve wave vector resolution along the chain 
axis. The background was separately measured at H = 
and subtracted from the data shown. In all scans the 
magnetic peak width was determined to be resolution- 
limited. Thus, to within the experimental wave vector 
resolution, the high-field phase in NDMAP for H \\ c is 
characterized by true 3-dimensional long-range order. 

To determine the high-field spin arrangement, inte- 
grated intensities of 17 magnetic Bragg reflections were 
measured at T = 25 mK and H — 5.8 T using a stan- 
dard diffraction configuration (no collimators). The ob- 
served intensity pattern is well reproduced by the simple 
model for the magnetic structure illustrated in the insert 
of Fig. |21 This spin arrangement is the same as previously 
seen for H||[01, l]i The collinear model is clearly over- 
simplified, and the actual structure should be a canted 
one, with all spins slightly tilted towards the field and 
thus producing a net magnetization along that direction. 
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FIG. 3: Measured peak intensity of the (0, 1.5, 0.5) magnetic 
Bragg reflection in NDMAP as a function of magnetic field 
applied along the crystallographic c axis (symbols) . The solid 
line is a power-law fit to the data, that should only be viewed 
as a guide for the eye. Inset: proposed model for the spin 
structure in the high-field phase. Only the staggered magne- 
tization probed by neutron diffraction is shown. In addition, a 
weak "ferromagnetic" tilt of all spins along the field direction 
is expected. 
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FIG. 4: Elastic scans across the (0, 1.5,0.5) magnetic Bragg 
reflection in NDMAP at T = 30 mK and H — 5.8 T show 
resolution-limited peaks in all three directions. Solid lines are 
Gaussian fits to the data. 



In our experiment we can detect only the staggered part 
L of the magnetic moment which lies, within the exper- 
imental accuracy, along the crystallographic a axis, per- 
pendicular to the field direction. The local spin direc- 
tions are opposite on sites related by (0, 1, 0) and (0, 0, 1) 
translations, and the same on sites related by a trans- 
lation along (1,0,0). Comparing the experimentally de- 
termined intensities of nuclear and magnetic reflections 
provides an estimate for the total staggered magnetiza- 
tion per site: L — 0.9(1) /ib- This experimental value 
is about half of the classical sublattice magnetization for 
an ordered S — 1 system. 



B. Inelastic scattering 

1. Constant-q data 

The field dependence of the gap energies was mea- 
sured in a series of energy scans performed at the ID 
AF zone-center (0, fc, 0.5). The scans corresponded to a 



fixed momentum transfer parallel to the chains qu = n/c. 
The wave vector transfer perpendicular to the chains 
q±_ = 2nk/b was varied in the course of the scan to satisfy 
geometrical restrictions imposed by the construction of 
the horizontal-field magnet. The background for all scans 
was measured at H — and H = 6 T, away from the 
ID AF zone-center, at q = (0,fc,0.33). Apart from the 
expected elastic contribution due to incoherent scatter- 
ing, the background was found to be energy-independent 
and about 1.5 counts/min. Typical background-removed 
scans collected at H = 1 T< H c , H = 4.1 Tw H c and 
H = 6 T> H c are shown in Fig. [S] (symbols). 

Scans collected at all fields are combined in the 3D 
plot shown in Fig. The data were analyzed using a 
single-mode cross section function similar to that used in 
Refill 

Sfau) « |/( g )| 2 p /~ C ° S(qC) x 

X [S(tlUJ — hujj^q) + S(hLO + fiUJj.q)] (1) 

(?k^ q ) 2 = A| +w 2 sin 2 (qc). (2) 
Here j = 1,2,3 labels each of the three excitation 
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FIG. 5: Typical background-subtracted constant-^ scans 
measured in NDMAP up to 1.75 meV energy transfer for dif- 
ferent values of magnetic field applied along the crystallo- 
graphic c axis (symbols). The two peaks are the lower-energy 
members of the Haldane excitation triplet. The third mode 
has a larger gap and is outside the shown scan range. Solid 
lines are fits to the data using a single-mode model cross sec- 
tion, as described in the text. 



branches, v is the spin wave velocity, f{q) is the mag- 
netic form factor for Ni 2+ , and Aj are the gap energies. 
The intensity prefactors Pj depend on both the matrix 
elements between the ground state and the single-mode 
excited states and on the polarization of the latter. The 
model cross section was numerically convoluted with the 
calculated spectrometer resolution function. The gap en- 
ergies and prefactors for each mode were refined using a 
least-squares routine to best-fit the data collected at each 
field. The spin wave velocity was fixed at v = 6.5 meV, 
as determined previously for H = 0.- The resulting fits 
are represented by the solid lines in Figs. |5] and H3 The 
shaded areas in Fig. [S] are partial the contributions of 
each mode. The obtained field dependence of the en- 
ergy gap is plotted in the (H, hu>) plane of the 3D plot 



in Fig. EJand, in more detail, in the top panel of Fig. 
We have observed only two lower-energy members of the 
Haldane excitation triplet; the third mode has a larger 
gap and is outside the shown scan range. The field de- 
pendence of the upper two triplet modes was observed in 
recent ESR measurements^ However, the lowest triplet 
mode was not observed in Ref. 4 since those measure- 
ments were done at much higher temperature T = 1.5 K 
where the lowest mode becomes strongly damped, the 
situation similar to that encountered in our early exper- 
iments in the H \\ a geometry which were also done at 
high temperature and failed to observe the reopening of 
the gap in the lowest mode. 



2. Constant-E scans 

As mentioned in the introduction, theory predicts that 
the spectrum at H > H c in the ideal AS geometry should 
lose its single-mode character. The sharp magnon excita- 
tions are expected to be replaced by a diffuse continuum 
of states with a lower bound following the magnon disper- 
sion curve. Even in this case, the continuum is singular at 
the lower bound and may be difficult to distinguish from 
a single-mode excitation smeared effects of experimental 
resolution. 

In search for any deviations from the single-mode pic- 
ture we performed constant-^ scans at 1.8 meV (Fig. 0) 
and 1.2 meV (not shown) at H = 2 T< H c , H = 4 T« H c 
and H = 6 T> H c . These particular energies were se- 
lected to avoid both gap energies at all three field values. 
A constant background was assumed for each scan. At 
each field, a single-mode profile was simulated using the 
model cross-section function described above, with pa- 
rameter values obtained in the analysis of const-<?|| data. 
We found that, to within experimental error, the single- 
mode model (solid lines in Fig. [7J reproduces all mea- 
sured const-i? scans very well, below, at and above H c . 
The observed variation of scan shape is due to changes in 
the gap energies of the two lower modes, shown as shaded 
areas in Fig. No features beyond those given by the 
single-mode approximation could be detected with the 
resolution of the present study. 



IV. ANISOTROPIC GAPPED QUASI-1D SPIN 
SYSTEM IN STRONG MAGNETIC FIELD: 
THEORY 

We now turn to developing a theoretical model of the 
high-field state. Our goal is a semi-quantitative effective 
field theoretical description, backed by a simple physi- 
cal picture, yet capable of consistently reproducing all 
the available experimental data on NDMAP. The latter 
implies that the model should account for both neutron 
and ESR measurements,"^ work both above and below 
the critical field, and apply in the case of arbitrary field 
orientation. 
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FIG. 6: Background-subtracted constant-^ scans measured in NDMAP for different values of magnetic field applied along the 
crystallographic c axis. Solid lines and shaded areas are fits to the data using a model cross section function, as described in 
the text. The black and white circles in the (H, hu) plane are the measured field dependence of the two lower energy gaps in 
NDMAP. The connecting solid lines are guides for the eye. 



A. A single anisotropic Haldane chain in a field 

1. Existing models 

In the early 90s, several phenomenological field- 
theoretical descriptions of the high-field regime in 
the anisotropic Haldane chain were proposedii§iiL2&i2& 
Affleck 17-25 proposed a theory based on coarse-graining 
the 0(3) nonlinear sigma model (NLSM). 1 ^ Technically 
coarse-graining leads to relaxing the unit vector con- 
straint of the NLSM, so that one has a theory of un- 
constrained real vector bosonic field <p. The </? 4 -type in- 
teraction was added to ensure stability in the high-field 
regime, and the anisotropy was introduced just by as- 
suming three different masses A a for the three field com- 
ponents, so that the resulting Lagrangian had the form 

£ = ^-{(d t <p + Hx<p) 2 -v 2 (d x <p) 2 

- ]>>^ 2 -A(<p 2 ) 2 }. (3) 

a 

For H > H c the ground state acquires a nonzero stag- 
gered magnetization L — (ip) and a uniform magneti- 
zation M oc (if x tp). This model captures the ba- 
sic physics involved, but is known to suffer from several 
drawbacks. Because of the too simplistic way of intro- 
ducing the anisotropy, the predicted values of the critical 



field (He = Aq, for the field if directed along one of the 
symmetry axes e a ) disagree with the results of pertur- 
bative treatmenlii 5 * 2 ^ 2 ^ as well as with the experimental 
data on the behavior of gaps as functions of the applied 
field in NENBSa and NDMAP^ 

TsvelikaS proposed a different theory which stems from 
the integrable Takhtajan-Babujian model of a S = 1 
chain and involves three Majorana fields with masses A a . 
The theory yields the critical field value H^ = ■^/A ( gA 7 
which coincides with the perturbative formulas f 15 - 27 - 28 . 
This model was rather successful for the description of 
the field dependencies of the gaps below H c in NENP, 29 
and is expected to yield a correct critical behavior at 
H — > H c . However, when the high- field neutron data 
on NDMAP in the H || a geometry became available^ 9 
it turned out that Tsvelik's theory, apart from overesti- 
mating the H c value considerably, predicts no change of 
slope for the two upper magnon modes at H = H c , in 
complete disagreement with the experimental data. One 
may conclude that though this model correctly describes 
the behavior of the low-energy degrees of freedom near 
the critical field, but fails to describe the behavior of 
high-energy modes above H c . 

Mitra and HalperinSS. have modified Affleck's bosonic 
Lagrangian in order to reproduce Tsvelik's results for the 
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FIG. 7: Background-subtracted constant-iJ scans measured 
in NDMAP for different values of magnetic field applied along 
the crystallographic c axis (color symbols). Solid lines and 
shaded areas are simulations using a single-mode cross section 
function and parameters obtained in the analysis of const-Q|| 
scans, as described in the text. 



gaps. It turns out that changing the first term in to 
i- £ {d^ a + ^2(A^/A a )^e afiy H^} 3 (4) 
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exactly reproduces the results of Ref. |l6| for the field de- 
pendencies of the gaps below H c . Above H Cl the pre- 
dicted field behavior of the gaps is different from that 
of Ref. El and is in a reasonable qualitative agreement 
with the experimental data on NDMAP in the H || a 
geometry^ However, apart from the fact that the rea- 
sons for the postulated modification remain unclear, 
the theory still has one fundamental flaw: it predicts 
that the staggered moment at H > H c is directed along 
the magnetic hard axis c for H \\ a and along the in- 
termediate axis b for H \\ c. This result is not only 
counter-intuitive but also contradicts to the diffraction 
experiments on NDMAP- which show that the ordered 
moment, in complete analogy to the classical picture for 



an ordered antiferromagnet, always lies in the easy plane 
along the most easy axis perpendicular to the field, i.e., 
along the b axis for H \\ a and along the a axis in the 
H II c case. 



2. An improved model 

We see that none of the previously known models pro- 
vides a consistent description of the experimental data. 
We use a different, more general approach, based on the 
model proposed in Ref. 30 for dimerized S =1/2 chains 
and S — 1/2 ladders, known to be in the same univer- 
sality class as S = 1 Haldane chains. This model was 
recently applied with great success to the description of 
the INS data for NDMAP in the H \\ a geometry^ and 
to ESR experiments in both geometries^ 

We first illustrate the general features of the theory 
on the example of the alternated S = i chain consisting 
of weakly coupled anisotropic dimers, described by the 
Hamiltonian 

W = ^2 JaS^n-iS^n + ^^iJ' (S2n- S2n+l)— H ■ S n }, (5) 
na n 

where < J' <C J. Throughout the rest of this sec- 
tion, it is implied that the magnetic field is measured in 
energy units, i.e., H i— > Q[IbH unless explicitly stated 
otherwise. For the derivation of the effective field theory 
it is convenient to use the dimer coherent states^ 

\A, B) = (1 - A 2 - B 2 )V2| S ) + £(A, + iBAfa), (6) 

3 

where the singlet state |s) and three triplet states \tj) are 
given bySi 

k> = -^(IU>-UT», = 7^(1 U> + lit)), 
IO = -^(1 TT) - I ii», I*,) = -J=(| TT) + 1 ii», 

and A, B are real vectors which are in a simple manner 
connected with the magnetization M = (Si + S2) and 
sublattice magnetization L = (Si — S2} of the spin dimer: 

M=2(lxB), L = 2(1-A 2 -B 2 ) 1/2 A. (7) 

The configuration space is the inner domain of the unit 
sphere A 2 + B 2 < 1 in i? 6 , with additional identification 
of the opposite points on the sphere, and the measure is 
defined as 67t _3 gL4 dB. 

We will assume that we are not too far above the criti- 
cal field, so that the magnitude of the triplet components 
is small, A, B <C 1. Then the effective Lagrangian den- 
sity in the continuum limit takes the following form: 

C = -2hB-d t A-^J'£ 2 (d x A) 2 -Y,{m l A 2 +m l B 2 } 



2H ■ (A x B) - U 4 (A,B), 



(8) 
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where mt 



4 \ e ijn 



m, — J', € is the lattice constant, and fhi 
J n ). The fourth-order term 



Vi(A, B) = \{A Z Y + X^A'B 2 ) + \ 2 {A ■ B) 2 , (9) 

where A = J', Ai = 2J', A2 = —J' in the present case. 
The spatial derivatives of B are omitted in (JHJ because 
they appear only in terms which are of the fourth order in 
A, B. Generally, we can assume that spatial derivatives 
are small (small wave vectors), but we shall not assume 
that the time derivatives (frequencies) are small since we 
are going to describe high-frequency modes as well. 

The vector B can be integrated out, and under the 
assumption A <C 1 it can be expressed through A as 
follows: 



B = QF, F = -hd t A +(H x A) 
n _ % \ AiAj 

^iij — — ^2 



vrii vnivrij 



(10) 



After substituting this expression back into © one ob- 
tains the effective Lagrangian depending on A only: 



rrii J 



(11) 



2-^(H x A)idtAi - U 2 (A) - Ui(A,d t A), 
rrii 



where Vi = U J'fhi/2 are the characteristic velocities in 
energy units, and the quadratic and quartic parts of the 
potential are given by 



U 2 {A) = miA 2 



-(H x A)\ 



(12) 



U 4 {A,d t A) = \(A 2 ) 2 + \ 1 A 



Ff + A 2 



A{Aj 
fhifhj 



-FiFi 



Note that the cubic in A term in (|10fl must be kept since 
it contributes to the U4 potential. 

Having in mind that the alternated S — \ chain and 
the Haldane chain belong to the same universality class, 
one may now try to apply this model in the form (|11I12|) 
to the Haldane chain system such as NDMAP, treating 
the velocities Vi and interaction constants m t , fhi, A, and 
Ai,2 as phenomenological parameters. 

One can show that the Lagrangian l|llfl contains theo- 
ries of AffleckiiS^ and Mitra and Halpcrm 26 as particular 
cases. Indeed, restricting the interaction to the simpli- 
fied form with A 1.2 = and assuming isotropic velocities 
Vi = v, one can see that the Affleck's Lagrangian © 
corresponds to the isotropic B-stiffness fhi — fh, while 
another choice m, = m, yields the modification 

For illustration, let us assume that H || z. Then the 
quadratic part of the potential takes the form 



U 2 = {m x 



H -)Al + (my - 5-)A$ + m z A 2 z , (13) 

Tfl x 



and the critical field is obviously H c = 
min{(m x rhy) 1 ^ 2 , (mythx) 1 / 2 }. At zero field the three 
triplet gaps are given by Aj = (mifhi) 1 ^ 2 . Below H c the 
energy gap for the mode polarized along the field stays 
constant, tko z = A z , while the gaps for the other two 
modes are given by 



1(^1 

2 V x 



A 2 ) + H 2 



± 



(A 2 X ~ A 2 ) 2 + H 2 (m x + m y )(m x + m y ) 



(14) 
1/2 



Below H c the mode energies do not depend on the inter- 
action constants A^. 

It is easy to see that in the special case nij = fhi, the 
above expression transforms into 



flOJZ 



\(A X 



Ay)± 



Ay) 2 +H' 



1/2 



(15) 



which exactly coincides with the formulas obtained in 
the approach of Tsvelikj^ and also with the perturba- 
tive formulas of Rcf. Il5ll27ll2l and with the results of 
modified bosonic theory of Mitra and Halperim 2 - The 
peculiarity of this special choice of parameters is that 
above H c = (uijUij) 1 ' 2 the most negative eigenvalue of 
the quadratic form l|13|) corresponds not to the compo- 
nent of A along the easy axis in the (xy) plane, as one 
would intuitively expect, but to the component along the 
harder axis. For instance, if x is the easy axis, m x < m v , 
then the most negative coefficient will be that at A y . 
This leads to the above-mentioned problem with counter- 
intuitive direction of the ordered moment in the theory 
of Mitra and Halperin. 

Generally, at H > H c one has to find the minimum of 
the static part of the potential and linearize the theory 
around the new static solution A = A^ . The equations 
for A^ have the form 

^2 ^Po A i + &{3i,mnAiA m A n = 0, (16) 
j imn 

where the matrices Q, A are defined as 

Qij = TTlidij — €ikn€jln-^—HkHi, 



kin 

a — r - 4- r 



J- ij,mn — A(5y S mn + (\16ij + A 2 ) Ci/tnE 



jlm- 



(17) 
HkHi 

rhifhj 



The magnon energies htu as functions of the field H and 
of the longitudinal (with respect to the chain direction) 
wave vector q can be found as three real roots of the 
secular equation 

det(M - (hu) 2 G - ihujC) = 0, (18) 

where symmetric matrices M and G are given by 



m 



y 



Mij = n tj + 



{Viqlf 
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ran 

Gij = l^--Xi(^r-) \s tj -X 2 -^-, (19) 

and the antisymmetric matrix C = R — R T is determined 
by 

^ = {i _Ai (t) }^ lHl 

A {0) r 1 1 i 

+ ^^^-^W^. (20) 
rrii ^ Tfij m n J 

Here in case of NDMAP the longitudinal wave vector q 
must be understood as counted from the ID Bragg point, 
ql i ^ 9||C— 7T. 

B. Inter-chain interactions 

Up to now, we have treated the problem as purely one- 
dimensional. In NDMAP, however, antiferromagnetic in- 
terchain interactions along the crystallographic b direc- 
tion lead to an observable transverse dispersion with the 
bandwidth of about 0.1 meVi^A This value is small com- 
pared to the magnon bandwidth along the chain axis 
(~ 7 meV), but constitutes approximately 20% of the 
lowest magnon gap, so that transverse interactions have 
to be taken into account if one aims at a quantitative 
description. 

It is straightforward to incorporate this effect into our 
formalism. Additional coupling of the form 

Jx$>*i-4w). (21) 

3 

where j labels chains along the transverse 6-direction, 
amounts to a renormalization of the A-stiffnesses m a i— > 
m a + 2J± cos(q±b). The minimum of magnon energies is 
reached at the 3D AF zone center, in our case at que = 
q\b = 7T. It is convenient to redefine the stiffness m a as 
its value at the zone center m a (q± — 7r/6), so that in the 
formulas l|lt)l) - (|20|l one just has to make the substitution 

m a I— > m a + 2Jj_(l + cos q±b). 

Secular equations (|18H then yield magnon energies for an 
arbitrary transverse wave vector transfer. 

C. Comparison with experiment 

The main advantage of the described model is that 
it can consistently reproduce all the experimental data 
currently available for NDMAP. At a first glance, it may 
appear over-parameterized, with 9 separate phenomeno- 
logical constants: A, Ax and A 2 . However, all 




H(T) 

FIG. 8: Symbols: measured field dependence of the energy 
gaps in NDMAP for a magnetic field applied along the crystal- 
lographic c-axis (top, this work) and a-axis (bottom, Ref.0). 
Lines: a global fit to the experimental data based on the 
Ginsburg-Landau type model outlined in Ref.|^and described 
in detail in the text. Note the different scales on the top and 
bottom figures. 

these parameters are relevant and can be almost uniquely 
determined from the measured field dependencies of the 
gap energies. Indeed, as discussed above, the indepen- 
dently measured zero-field gaps A a fix three relations 
= m a m a . The value of the critical field Hc^ for 
a field applied along the principal anisotropy axis a 
determines another three relations between parameters, 
namely (Hc^) 2 — min(m0m 7 , m 7 TOg). One can show 
that of those six equations for six stiffness constants m a 
and rh a only five are independent, so that all stiffness 
constants can be expressed through one of them (we have 
chosen m x for this role) . Finally, the interaction parame- 
ters A, Ai ; 2 control the behavior of the gap energies above 
the critical fields. In fact, the gaps depend only on the 
relative interaction strengths (Ai/A), (A2/A), so that the 
scale of A does not influence the expressions for the gaps 
and can be set deliberately (we have put A = 1). Thus, 
the knowledge of the H = gaps and of critical fields 
helps to fix five parameters, and one parameter turns 
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out to be irrelevant, so that one is left with only three 
parameters to fit the A(i?) curves. 

In analyzing the measured field dependencies of the 
gap energies in NDMAP, one has to keep in mind that 
both for the experiments described here, and for those 
reported in Ref. for H || a, the transverse wave vec- 
tor q± = 2itk/b is not constant, but varies as a func- 
tion of energy transfer hto. For H || c, q± is dic- 
tated by the geometry of the horizontal-field magnet, 
g ± 6/(27r) = 0.42 mcV -1 x (hu>). In the H \\ a experi- 
ment of Ref. qj_ was chosen to optimize wave vector 
resolution along the chains when using a horizontally fo- 
cusing analyzer, q±b/ (2tt) — 1.3 + 0.24 meV -1 x (hcu). In 
both cases the presence of inter-chain interactions must 
be taken into account explicitly, as discussed in the pre- 
vious subsection. The value of 4 J± was chosen to match 
the transverse dispersion bandwidth of 0.1 meV observed 
in NDMAP 24 at H = 0: Jj_ = 0.025 meV. 

A global fit to the data of both experiments is shown 
in Fig. [S] The coordinate axes were chosen along the 
local anisotropy axes for each Ni 2+ ion: y is parallel 
to the b axis, z is in the (a, c) plane and forms an 
angle of 16° with the c axis and x completes the or- 
thogonal set. The final parameters obtained in the fit 
are: m x = 0.50, m y = 0.71, m z = 4.76, rh x = 0.35, 
m y = 0.38, m z = 0.71, Ai/A = 0.17 and A 2 /A = -0.17. 
The solid lines in Fig.[S]never quite reach zero at H c , since 
they actually represent excitation energies at qy = n and 
transverse wave vector transfers matching those probed 
in the corresponding experiments. At the critical field 
these energies are non-zero due to non-zero dispersion 
perpendicular to the chain axis (see above discussion). 

Note that the fitted critical fields He = 5.5 T and 
= 3.4 T are some 15% smaller than observed in 
our neutron scattering experiments. Incidentally, these 
critical fields are in excellent agreement with those found 
in heat capacity 5 ESR measurements A The parameter 
values obtained in the present neutron study also agree 
nicely with those determined in the analysis of ESR res- 
onance frequencies. Finally, we have verified and would 
like to stress that the obtained parameter values yield 
correct directions of the ordered staggered moment at 
H > H c for both H \\ c and H \\ a experimental geome- 
tries: namely, the order parameter A is directed along the 
b and a axes in the H \\ a and H \\ c cases, respectively. 

V. CONCLUDING REMARKS 

Even for H \\ c NDMAP shows no signs of Luttinger 
spin liquid behavior, such as incommensurate correla- 
tions or breakdown of the single-particle spectrum. On 
the contrary, the system was shown to be antiferromag- 



netically ordered (a "spin solid" state) at high fields with 
an appreciable sublattice magnetization. This is accom- 
panied by a re-opening of the gap at H > H c . Overall, 
the observed field dependence of the excitation spectrum 
is qualitatively similar to that previously seen for H _L c& 
The reasons for the Luttinger spin liquid regime being 
unobservable are quite clear: the idealized axially sym- 
metric geometry can not be realized in NDMAP. This, 
to our opinion, may be mainly attributed to the effects 
of the 16° canting of the principal anisotropy axes rel- 
ative to the field direction, which define the strongest 
explicit breaking of the axial symmetry, although in- 
plane anisotropy and inter-chain interactions play a sig- 
nificant role as well. Perhaps in very strong external 
fields, when the Zeeman energy becomes large compared 
to any anisotropy and 3D effects, certain features of the 
Luttinger spin liquid may become accessible. In partic- 
ular, it has been recently argued^ that incommensurate 
correlation should emerge even in the axially asymmet- 
ric geometry, above a certain second critical field H*. 
Whether or not such an experiment is technically feasi- 
ble for NDMAP is currently unclear. 

In summary, the main impact of the present H \\ c mea- 
surements is to provide crucial quantitative data needed 
to evaluate the veracity of the various field-theoretical de- 
scriptions of magnetized anisotropic Haldane spin chains. 
Our conclusion is that of several existing theories, only 
the newly proposed model is robust enough to reproduce 
all the features observed experimentally for different field 
orientations. 
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